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ABSTRACT 

Aims. We use a combination of CFHTLS deep optical data, WIRCam Deep Survey (WIRDS) near-infrared data and XMM-Newton 
survey data to identify z > 1.1 clusters in the CFHTLS Dl and D4 fields. Counterparts to such clusters can not be identified without 
deep near-infrared data and as such the total of » Ideg^ of /, H and K, band imaging provided by WIRDS is an indispensable tool in 



o 
u 

(~ l , such work. 

^ |i Methods. Using public XMM X-ray data, we identify extended X-ray sources in the two fields. The resulting catalogue of extended 

I ' X-ray sources was then analyzed for optical/near-infrared counterparts, using a red-sequence algorithm applied to the deep optical 

O ' and near-infrared data. Redshifts of candidate groups and clusters were estimated using the median photometric redshifts of detected 

counterparts and where available these were combined with spectroscopic data (from VVDS Deep and Ultra-Deep and using AAT 

^ . AAOmega data). Additionally, we surveyed X-ray point sources for potential group systems at the limit of our detection range in the 

^, ' X-ray data. A catalogue of z > 1.1 cluster candidates in the two fields has been compiled and cluster masses, radii and temperatures 

have been estimated using the scaling relations. 

*^ , Results. The catalogue of group and cluster candidates consists of 15 z > 1.1 objects. We find several massive clusters (M > lO'^^M©) 

►^ . and a number of lower mass clusters/groups. Three of the detections are previously published extended X-ray sources. Of note is 

VQ ' JKSC 041 (previously detected via Chandra X-ray data and reported as a z = 1.9 cluster based on UKIDSS infrared imaging) for 

f^ ' which we identify a number of structures at redshifts of z = 0.8, z = 0.96, z = 1.13andz= 1.49 (but see no evidence of a structure at 

fvj , z = 1.9). We also make an independent detection of the massive cluster, XMMXCS J2215. 9-1738, for which we estimate a redshift 

w~, ■ of z = 1.37 (compared to the spectroscopically confirmed redshift of z = 1.45). We use the z ^ 1.1 catalogue to compare the cluster 

• ' number counts in these fields with models based on WMAP 7-year cosmology and find that the models slightly over-predict the 

C^^ ' observations, whilst at z > 1.5 we do not detect any clusters. We note that cluster number counts at z <; 1.1 are highly sensitive to the 

^^ , cosmological model, however a significant reduction in present statistical (due to available survey area) and systematic (due to cluster 

^^ . scaling relations) uncertainties is required in order to confidently constrain cosmological parameters using cluster number counts at 

high redshift. 

Key words, methods: data analysis surveys galaxies: clusters: general cosmology: large-scale structure of Universe 
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_C3 1. Introduction iKomatsu et"aLll2009l; IViMilinin et al.ll2009b . Recently, attention 

of theoretical cosmologists has been focused on the possibility 
Development of large-scale structure in the Universe is a sensi- ^f constraining different models of inflation through their pre- 
tive cosmological tool, which indepedently confirms the success dictions of departur es from Gaussianity in the primordial cur- 
of the ACDM cosmological paradigm of inflationary theory (e.g. y^t^re fluctuations (iDesiacaues et al.ll2009h . In addition, some 
coupling between dark matter and dark energy fluids are pos- 



* Based on observations obtained with MegaPrime/MegaCam, a joint sible, whi ch can give testa ble changes in growth of large scale 

project of CFHT and CEA/DAPNIA, at the Canada-France-Hawaii structure JBaldi et al.ll2010h . 

Telescope (CFHT) which is operated by the National Research Council /-^i^ ri- i.-u -j-^uj 

/T^Tn^x rA J \i. T .■» 5xT .■ 1 J o ■ J I'TT ■ jr.u Clusters of galaxies, which reside in the deepest gravita- 

(NRC) ol Canada, the Institut National des Science de 1 Univers ot the . , • , ,, ,, , • , , 

Centre National de la Recherche Scientifique (CNRS) of France, and Clonal potential wells, present an excellent cosmological probe 

the University of Hawaii. This work is based in part on data products as a tracer of the the peaks of structure formation. For example, 

produced at TERAPIX and the Canadian Astronomy Data Centre as Allen et al. (2008) used the cluster X-ray mass fraction (i.e. the 

part of the Canada-France-Hawaii Telescope Legacy Survey, a collabo- ratio of the baryonic mass density to the total mass density) to 

rative project of NRC and CNRS. place constraints on the cosmological mass density, Q,„, whilst 
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also providing an inde pendent constraint on t he cosmic accel- 
eration. Further to this, IVikhhnin et alJ (l2009h used the cluster 
mass function to measure structure growth to constrain the evo- 
lution of cosmic acceleration, measuring the cosmic acceleration 
parameter, w p, with a result o f wp = -1.14 + 0.21. 

Recently. i Jee et al.) ( l2009l) presented a w eak lensing analysi s 
ofthez = 1 .4 cluster, XMMU J2235. 3-2557 ('Mullis et al."2005'), 
measuring a cluster mass of x 8.5 x 10'''Mo. Based on WMAP 
5yr cosmology they predict the expected number of clusters of 
this mass at this redshift to be lO"-' in their lldeg^ survey. It is 
suggested by the authors that the discovery of such a rare cluster 
within the given survey area is indicative of either a higher value 
of cTg, or stru cture growth from non-ga ussian primordial density 
fluctuations. iJimenez & Verdd (l2009h investigate the effect of 
non-gaussianity on the number density of such massive clusters 
at the redshift of XMMU J2235. 3-2557 and conclude that mass 
determination of high-redshift, massive clusters can offer a com- 
plementary probe of primordial non-gaussianity. Further to this, 
the detection of the massive X-ray selected cluster, XMMXCS 
J2215.9-1738 ('Stanford_eLaD|2Q06l) at a redshift of z = 1.45 has 
been treated as another supportive argument to these findings 
(IJimenez & Verdell2009h . 

Uncertainty in the current measurements of erg is also 
hinted at from studies of the thermal Sunyaev-Zel'dovich 
(tSZ) effect from clusters. A number of papers have claimed 
lower than expected tSZ signals from clusters in Cosmic 
Microwave Back ground (CMB) maps (e.g . JLieu et al. 2006; 
Bielby & Sha nksl 120071; [Luekeretal.ll2O10t lOiego & Partridge, 

iarticular|D 



2010; Komatsu et al. 2010) and in parti cular lLueker et all (1201 Ol) 



claim this to be suggestive of a lower value of erg. From their 
measurements using the South Pole Telescope (SPT), they give 
constraints of erg, = 0.773 + 0.025 compared to WMAP 5-year 
constraints of erg = 0.812 + 0.026 from the WMAP-i-BAO-i-SN 
combination (although from WMAP alone a more consistent 
value of cTg = 0.796+0.036 is calculated). Although these differ- 
ences are only at the ~ Icr level, they can have significant effects 
when used in simulations of structure formation (in particular on 
the cluster number counts). From both tSZ and cluster number 
studies, it is evident that significant uncertainty still remains in 
the determination of erg, whilst these uncertainties can inhibit 
our ability to constrain non-gaussianity in the primordial density 
fluctuations. 

Critical to applying cluster observations to constraining the 
cosmological framework are the X-ray scaling laws that may 
be used to relate cluster properties such as mass, temperature 
and X-ray luminosity. From a simple model of pure gravita- 
tional collapse (referred to as the self-similar model) for ex- 
ample, power-law relations are predicted for the Lx - T and 
Lx - M cluster relations ( Kaiser 1986). Much work ha s gone 
into characterising these relations (e. g. Markevitch et al .1ll998t 
Reiprich & Bohringei^l2002l: IVikhlinin et al.„2006c .Rvkoff^ et al.1 
.2008; Comerford et al. 201Ci), with calibrations having been per- 
formed up to redshifts of z ~ I (Leauthaud et al. 2010). 

At the present time there are three promising ways to find 
a distant galaxy cluster - via the Sunyaev-Zel'Dovich Effect, 
through its extended X-ray emission, using the red-sequence 
cluster method with near-infrared data. Optical surveys are lim- 
ited to z ~ 1.1 as the 4000A break moves into the infra-red 
and can no longer b e used as an ident ification tool in optical 
bands. For example, lOlsen et al] (l2007l) presented a catalogue 
of galaxy clusters identified in all four CFHTLS Deep fields 
purely from the CFHTLS optical u'griz data (i.e. with no prior 
X-ray detections). The catalogue provides 169 cluster candidate 
detections with a maximum redshift of z = 1.2. Similarly, the 



Bon n Lensing, Optical a nd X-ray (BLOX) galaxy cluster sur- 
vey (iDietrich et al.ll20()/ ) provides a catalogue of weak-lensing 
and X-ray selected clusters in a number of fields covering a total 
area of 6.4deg^. In this case the limiting redshift of the catalogue 
is z ~ 1 .0 as their imaging data is limited to the BVRI optical 
bands, leaving a number of cluster detections in their catalogue 
with no available redshift estimate. The weak-lensing/X-ray ob- 
servations can not by themselves separate out the high-redshift 
cluster candidates and t heir spectros copic follow-up is a tedious 
(though rewarding e.g. Mullis et al.ll20 04) task. 

The Palomar Distant Cluster Survey (PDCS: IPostman et al.l 
1996) performed an optical selection of galaxy clusters over an 
area of 5.1deg''. This catalogue covered a range in redshifts of 
0.2 < z < 1.2 identifying 107 cluster candidates, but with no X- 
ray data provided no evidence of the potential wells tha t these 
potential clusters resided in. More recently, Finogueno v et al.l 
(2007) presented cluster detections across the Cosmic Evolution 
Survey (COSMOS) field. Based on initial X-ray detections 
from XMM data, they identify 72 clusters across an area of 
2deg" using optical + NIR data to determine photometric red- 
shift estimates. Their fu ll sam ple reaches a redshift limit of 
z ~ 1.3.|E isenhardt et alj (l2008h used the IRAC Shallow Survey 
(lEisenhardt et al.ii2004) to identify 106 z > 1 galaxy cluster can- 
didates over an area of 7.25deg^. Subsequent spectroscopic ob- 
servations of cluster members confirmed 12 of these candidates 
with further observations to follow. They determine initial clus- 
ter redshift estimates using photometric redshifts based on deep 
optical data from the NOAO Deep Wide-Field Surv ey (NDWFS; 
Jannuzietal. 2004), NIR data from FLAMEX dElston et al.1 
2006) and IRAC Shallow Survey [3.6] and [4.5] imaging. With 
this combination they successfully identify clusters up to red- 
shifts of z ~ 1 .5 , with the success o f the identification of a 
z = 1.41 cluster (I Stanford erani2005h . at the time the highest 
redshift spectroscopically confirmed galaxy cluster At z > 1.5 
however, they found their survey limited by the depth constraints 
of the optical data available, despite the capabihty of their IR 
data. This was followed by the discovery of the z = 1 .45 clus- 
ter XMMXCS J2215.9-1738 (located in the CFHTLS D4/LBQS 
2212-1759 field) by Stanford et al. (2006), which was initially 
detected in the X-ray with XMM data, with / and Ks band imag- 
ing an d with Keck DE IMOS spectroscopic follow-up. More re- 
cently, [Andreonetal] (|2009) presented a cluster candidate at an 
estimated redshift of z - 1.9 for an extended X-ray source, 
based on red-sequence analysis using deep NIR data from the 
UKIRT Deep Infrared Sky Survey. However, this remains un- 
confirmed by spectroscopic observations at the present time. The 
highest confirmed redshift cluster presently stands as the z = 
1.62 in the Subaru/XMM- Newton De e p Fiel d (S XDF) discov- 
ered in dependently by Papovich et al.l (1201 Ol) and iTanaka et al.l 
(l2010h . Both present red-sequence analysis and spectroscopic 
follow-up showing that the red-sequence at bright magnitudes 
is in place at least to z = 1.6. Recently the Spitzer Adaptation 
of the Red-Sequence Cluster Survey (SpARCS Muzzi n et al] 
2009; Wilson et al. 2009; Demai-co et al. 2010) has used the red- 
sequence method to identify z > 1 clusters in six fields covering 
a total area of ^ 50deg^. A number of these have now been suc- 
cessfully confirmed with spectroscopic follow-up. 

Deep NIR surveys are required to identify clusters at z > 1 , 
but are only available on a limited area of the sky. A number 
of such surveys have X-ray coverage s uflicient for high-redshift 
cluster search ("Finogu enov et al.l 1201 0) and in this paper we 
present an overview of a promising combination between one 
of the best deep NIR data-sets avai lable to date, the WIRCam 
Deep Survey (WIRDS: lBielbv et al.ll!n Prep) , with the CFHTLS 
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Deep optical data and X-ray data provided by XMM-Newton. 
Subsequent papers are planned, which will present the low- 
redshift detections and look at the properties of cluster members 
in the cluster sample. In section|2] we describe the data-sets used 
in this survey. Section [3] describes the methods used for cluster 
detection and redshift estimation, section |4] presents our high-z 
(z > 1.1) cluster candidates and section |5] presents the analysis 
of the high-z catalogue data-set. We conclude in section|6] 

Unless otherwise stated we use a flat ACDM cosmology with 
fi,„ - 0.25 and //q = 72 kms^'Mpc"', whilst magnitudes are 
given in the AB system. 



abundance and a redshift of 0.2. We note that in reconstructing 
the properties of the identified groups and clusters of galaxies as 
well as in the modelling of survey z-dependent characteristics, 
we implement the exact corrections, based on the source spectral 
shape (as defined by the expected temperature of the emission) 
and the measured redshift of the system. 

After the background has been estimated for each observa- 
tion and each instrument separately, we produce the final mosaic 
of cleaned images and correct it for the mosaic of the exposure 
maps in which we account for diff'erences in sensitivity between 
pn & MOS detectors. 



2. Data 

2.1. X-ray Data 

A description of t he XMM-Newton observatory is given by 
iJansen et al.l (l200lh . In this paper we use the data collected 
by the Euro pean Photon Imagin g Cameras (EPIC): the pn- 
CCD cam era 'Stru der etal] ( |2001|) and the MOS-CCD cameras 
dTurner et al. 200 15~^ 

XMM-Newton observations of the CFHTLS Deep fields are 
goals of different surveys. The Dl field has been covered within 
the XMM-LSS survey (PI M. Pierre). The D2 field is located 
within the COSMOS field and has the deepest X-ray coverage to 
date (PI G. Hasinger). The D3 field is also one of DEEP2 fields 
and the X-ray data has been taken there as a part of AEGIS sur- 
vey. The XMM observations of D4 field were obtained through 
a proposal of the AEGIS team (PI K. Nandra). 

To date most of these data have been publicly available 
through the ESA XMM-Newton Science Archive (XSAQ) and 
we have processed the X-ray data for all CFHTLS Deep fields. 
The coverage of the D3 by XMM is the smaflest of all (0.2deg^), 
since most AEGIS data has been obtained using Chandra and 
that program is still ongoing. The results of our survey for 
COSMOS has already been use d for t he identificatio n of X- 
ray clus ters in iFinoguenov et al.l (l2007h and iFinoguen ov et al. 



(In Prep ), as well as c luster science exploration flGiodini et al 



[2009; LeauthaudetalJbOlO) . As the COSMOS field is treated 
in depth in these other papers and the coverage of the D3 field is 
comparatively small, we focus here on the Dl and D4 fields. 

The XMM coverage of both D 1 and D4 are comparable with 
the entire area of these two fields being observed to 10-20 ksec 
depth, thus presenting a homogeneous dataset. All X-ray obser- 
vations that we have used in this work are listed in Table [1] 

The initial data processing has been done using the 
XMMSAS version 7.1 (Watson etal. 2001; Kirsch et al. 2004; 
fSaxton et al.l l2005). Upon creating the calibrated event files, we 
perform a more conservative removal of time inter vals affected 
by sola r flares, following the procedure described in lZhang et al.l 
(l2004l) . In order to increase our capability of detecting ex- 
tended, low surface brightness featu res, we have app l ied th e 
four-stage background subtraction of IFinoguenov et all (l2007l) . 
We also check for a high ba ckground that can be present in the 
MOS chips ( Kunt z & Snow den 2008), identifying and remov- 
ing from further analysis the following hot chips: MOSl chip 
4 in ODE 0404960101 and MOS2 chip 5 in ODFs 021049001, 
0404960101, 0404960501. The resulting countrate-to-flux con- 
version in the 0.5-2 keV band excluding the lines is 1 .59 x 10 '^ 

for pn and 5.41 x 10"'^ for each MOS detector, calculated for the 

■ 1 

source spectrum, corresponding to the APEC (Smith et al. 2001) 

model for a collisional plasma of 2 keV temperature, 1/3 solar 



2.2. Spectroscopic Redshifts 

In the CFHTLS Dl field, we use spectro scopic redshifts 
from the VVDS Deep dLe Fevre et alj l2005h and Ultra-Deep 
(ICucciati et al.lln Preph spectroscopic samples. The VVDS Deep 
sample is available publicly and consists of 898 1 spectroscopi- 
cally observed objects over an area of 1.2deg^ in the CFHTLS 
D 1 field. It consists of a magnitude limited sample with a limit of 
/ < 24 and samples a redshift range of < z < 5. The Ultra-Deep 
sample consists of ~ 1500 spectra over an area of a 0.14deg 
and covers a magnitude range of 22.5 < Iab < 24.75. Both of 
the VVDS spectroscopic catalogues attribute each object a flag 
based on the identification. These range from 1 to 4 with 1 being 
most unreliable and 4 being most reliable. In addition a flag 9 is 
given to objects identified based on a single emission line. 

Spectroscopic redshift data in the D4 field were ob- 
tained using the AAOmega instrument at the Anglo- Australian 
Observatory (AAO) as part of a program to provide optical spec- 
troscopy of X-ray point-sources in the CFHTLS ( StaUn et aP 
2010). Observations were taken on the dates 25th - 27th 
September 2006 (program ID: 06B/027, PI: C. S. Stalin) and 
11th - 13th September 2007 (program ID: 07B/026, PI: P 
Petitjean) with central coordinates of 22:15:30 and -17:44:00. 
The field of view of the AAT with the AAOmega instrument is 
defined by a diameter of 2° . The observations were made with in- 
dividual exposure times of 1800s during the 2006 run and 1680s 
in the 2007 run. In order to maximize the wavelength cover- 
age observations were made using two different set-ups, the first 
with the 580V grism (with a central wavelength of 4790A) and 
the second with the 385R grism (with a central wavelength of 
7200A). Both grisms have a MOS resolution of i? = 1300 with 
a dispersion of 1.033 A/pixel and 1.568A/pixel for the 580V and 
385R grisms respectively. 

Reductions were made using the AAO 2DFDR software 
to perform flat-fielding and arc-lamp wavelength calibration. 
Redshift identification was then performed on all objects indi- 
vidually by eye and using galaxy templates to determine galaxy 
redshifts. In total from the CFHTLS D4 AAOmega specti'o- 
scopic data, we have redshifts for 492 objects. 139 of these are 
QSOs, 52 are stars and 301 are galaxies, all at magnitudes of 
/ < 22.5. 



2.3. CFHTLS Deep Optical Data 

All optical data used in this work was taken on the Canada- 
France-Hawaii Telescope (CFHT) as part of the CFHT Legacy 
Survey (CFHTLS). Specifically, we utilize the T0006 release of 
the Deep survey in the CFHT Dl and D4 fieldfl The CFHTLS 



http://xmm.esac.esa.int/xsa/ 



http://terapix.iap.fr/cplt/T0006-doc.pdf 
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Table 1. X-ray Observations. 



Name 


Obs. ID 


R.A. 


Dec. 


PI 


Filter" 




Net Exp (s) 








(J2000) 






pn 


MOSl 


MOS2 


XMM LSS^ 


0037980201 


36.00000 


-3.83333 


M. Pierre 


Thin I 


7780.1 


11825.6 


12733.5 


XMM_LSSJ 


0037980301 


36.33334 


-3.83333 


M. Pierre 


Thin I 


7762.5 


I3I26.6 


13128.8 


XMM_LSS.4 


0037980401 


36.66666 


-3.83333 


M. Pierre 


Thin I 


3029.7 


5858.6 


6011.8 


XMM LSS_5 


0037980501 


37.00000 


-3.83333 


M. Pierre 


Thin I 


10294.4 


14390.3 


14906.3 


XMDSOM_l 


0109520101 


35.83333 


-4.16667 


K. Mason 


Thin I 


I857I.5 


24756.3 


24768.6 


XMDSOM_2 


0109520201 


36.83334 


-4.83333 


K. Mason 


Thin I 


17389.9 


24087.9 


24150.9 


XMDSOM_3 


0109520301 


36.50000 


-4.83333 


K. Mason 


Thin I 


15097.4 


20780.1 


21299.1 


XMDSOM_5 


0109520501 


35.83333 


-4.83333 


K. Mason 


Thin I 


17300.5 


23176.2 


23186.8 


XMDSOM_6 


0109520601 


35.66667 


-4.50000 


K. Mason 


Thin I 


15336.5 


21098.1 


21768.6 


XMDSSSC.l 


0111110101 


37.00000 


-5.16667 


M. Watson 


Thin I 


1 2051.0 


20382.5 


19159.7 


XMDSSSC^ 


0111110201 


36.66666 


-5.16667 


M. Watson 


Thin I 


2142.7 


5907.8 


7307.9 


XMDSSSCJ 


0111110301 


36.33334 


-5.16667 


M. Watson 


Thin I 


16788.3 


21222.0 


20939.6 


XMDSSSC.4 


0111110401 


36.00000 


-5.16667 


M. Watson 


Thin I 


I832I.3 


26235.8 


26616.6 


XMDSSSC^ 


0111110701 


36.66666 


-5.16667 


M. Watson 


Thin I 


5857.6 


11226.3 


10820.7 


MLS-1 


0112680101 


36.83334 


-4.16667 


M. Turner 


Thin I 


20412.9 


23226.3 


23846.5 


MLS-2 


0112680201 


36.50000 


-4.16667 


M. Turner 


Thin I 


6390.5 


8159.4 


8148.5 


MLS-3 


0112680301 


36.16667 


-4.16667 


M. Turner 


Thin I 


I684I.3 


21034.0 


20375.4 


MLS-5 


0112680401 


37.00000 


-4.50000 


M. Turner 


Thin I 


18032.3 


22559.3 


22823.7 


MLS-8 


0112680501 


36.00000 


-4.50000 


M. Turner 


Thin I 


I4I00.4 


20720.2 


20728.6 


MLS-7 


0112681001 


36.33334 


-4.50000 


M. Turner 


Thin I 


I8I99.3 


20108.7 


20330.6 


MLS-6 


0112681301 


36.66666 


-4.50000 


M. Turner 


Thin I 


9589.7 


15886.8 


15903.5 


XLSSJ022404.0-04132 


0210490101 


36.00834 


-4.18861 


L. Jones 


Medium 


61458.3 


73743.5 


71108.0 


XMM-LSS.4 


0404960101 


36.66666 


-3.83333 


M. Pierre 


Thin I 


3752.0 


15728.5 


15341.6 


XMDSOM_4 


0404960501 


36.16667 


-4.83333 


M. Pierre 


Thin I 


7480.8 


9590.4 


10312.9 


LBQS2212-1759 


0106660101 


333.88196 


-17.73492 


J. Clavel 


Thinl 


48139.6 


53831.9 


54507.2 


LBQS2212-1759 


0106660201 


333.88196 


-17.73492 


J. Clavel 


Thinl 


30559.4 


44450.2 


44544.6 


LBQS2212-1759 


0106660501 


333.88196 


-17.73492 


J. Clavel 


Thinl 


4391.1 


7675.9 


7775.8 


LBQS2212-1759 


0106660601 


333.88196 


-17.73492 


J. Clavel 


Thinl 


69304.4 


85376.6 


86686.9 


CFHTLSJD4_1 


0505460101 


334.22229 


-17.92925 


K. Nandra 


Medium 


15579.8 


24137.6 


24067.8 


CFHTLSJD4_2 


0505460201 


333.49088 


-17.87772 


K. Nandra 


Medium 


16975.0 


24725.6 


22828.1 


CFHTLSJD4_3 


0505460301 


334.14038 


-17.45453 


K. Nandra 


Medium 


21628.4 


27582.6 


27287.2 


CFHTLSJD4_4 


0505460401 


333.64246 


-17.44842 


K. Nandra 


Medium 


4293.1 


13787.9 


18142.2 


CFHTLSJD4_5 


0505460501 


333.85226 


-18.09394 


K. Nandra 


Medium 


11959.0 


19149.9 


19352.9 



Note that in all observations the same filters were used for the pn, MOSl and M0S2. 



Table 2. Photometric depths in the CFHTLS optical and WIRDS 
NIR imaging. 



J 



H 



DI 
D4 



27.0 
26.5 



26.8 
26.7 



r I z .. . 

(50% point-source completeness) 

:^7-^ ^r^-T^ ^r^^ 24J 247" 



K. 



26.3 
26.3 



25.3 
25.1 



25.5 
25.4 



24.7 



25.1 24.6 24.6 



Deep data consists of imaging with u*griz filters. 50% point- 
source completeness limits for this data are given in table|2] 



2.4. WIRDS Near-Infrared Data 

In addition to the optical data from the CFHTLS, we have used 
near-infrared J, H and Ks dee p imaging from the WIRCam 
InfraRed Deep Survey (WIRDS. JBTelbv et al-llln^epl) . WIRDS 
gives coverages of ~ 0.6deg and ~ 0.4deg in the CFHTLS D 1 
and D4 fields respectively. Again, 50% point-source complete- 
ness limits are given in table |2] 

We construct matched catalogues for both fields incorporat- 
ing all eight available imaging bands. The WIRDS and CFHTLS 
images are all provided with the same pixel-scale of 0.186" and 
the same area in each field, allowing catalogues to be simply ex- 
tracted using SExtractor in dual-image mode. We therefore pro- 
duce two catalogues for each field, the first using a gri x^ image 



as the detection image and the second using the Ks image as 
the detection image. The gri x^ images were in turn produced 
using SWARP in x" combination mode using the CFHTLS g, 
r and i band images in each field. The x^ mode produces an 
optimized combination of the input images where the output is 
effectively the probability of a given pixel being part of the sky- 
distributi on, based on a re duced x^ with the background distri- 
bution (SzalaxeFalIT99^. We note that all images have approx- 
imately the same seeing quality of ~ 0.6" - 0.7" and as such we 
do not perform any smoothing of the images to match FWHMs 
before running SExtractor. Photometric calibration of the in- 
frared observations was performed by matching to the 2MASS 
photometric data in these fields. 



2.5. Photometric Redshifts 

We deriv ed photometric redshifts using the Le Phare code 
dArnouts et al. 2002; Ilbert et al. 2006) with nx^ template-fitting 
method appHed to the CFHTL SAVIRDS 8 - band photometry 
(u*grizJHKs). Templates from .PoUetta et al] (l2007h were used 
in combination w ith the additional blue-galaxy templates of 
Ilbe rt et aP (1200?). The photo-z were estimated using the me- 
dian of the probability distribution function (PDFz) rather than 
the minimum of the x^ distribution. We show the comparison 
of the photometric redshifts to available spectroscopic redshifts 
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Fig. 1. Comparison of the spectroscopic and photometric red- 
shifts in the CFHTLS Dl and D4 fields. All photometric red- 
shifts are taken form our 8-band WIRDS photometric redshift 
catalogues. Spectroscopic redshifts are from the WDS Ultra- 
Deep (23 < / < 24.75) sample for the Dl data (black squares) 
and the AAOmega spectra described in s ection 12.21 (open blue 
circles) and 5 spectroscopic redshifts from lStanford et al] (l2006l) 
(filled red circles) for the D4 data. 



in the two fields (see section 12721 for details of the spectroscopic 
data) in Fig.[T] The black squares show the comparison between 
our photometric redshifts and the WDS Ultra-Deep data in the 
CFHTLS Dl field, whilst the blue circles show the comparison 
with AAOmega spectroscopic redshifts in the D4 field and the 
red filled circles s how the results for galaxies identified in the 
z - 1.45 cluster of [Stanford et al.l (12006) . We find accuracies of 
the photometric redshifts of cr az/(i+z) - 0.032 with a failure rate 
of 2.5% when comparing the photometric redshifts to the WDS 
Deep data in the Dl field. Comparing with the WDS Ultra- 
Deep spectroscopic redshifts, we find o'aj:/(i+,) = 0.037 with a 
failure rate of 1%. These figures are derived across the redshift 
range < z < 2 and using only flag 3 and 4 objects from the 
spectroscopic datasets. We note that we see a slight systematic 
under-prediction of redshifts at z > 1 . As discussed, the PSF in 
the optical and infrared images are relatively consistent and we 
have not performed any PSF matching. However, the minor dif- 
ferences in the PSF are the likely cause of the observed system- 
atic offset. Although it will cause us to marginally underestimate 
cluster redshifts (by Az ~ 0.05), it should not affect our ability 
to identify galaxy over-densities, whilst the cluster redshifts may 
be more precisely identified using spectroscopic follow-up. 



3. Cluster Detection 

3.1. Extended X-ray Source Detections 



We use the prescription of lFinoguenov et all (1201 Ol) for extended 
source detection, which consists of removing the PSF model for 
each detected point source from the data before applying the ex- 
tended source search algorithm. We in clude an additiona l step i n 
background refinement performed bv iFinoguenov et aP (1201 Ol) . 
in which we repeat the background estimate steps for each detec- 
tor and pointing. This refines the definition of the area used for 
the background evaluation, which is obtained through the analy- 
sis of the final mosaics. In particular, the area free from expected 
contamination of sources is refined based on the PSF model im- 
age, in which we excise zones where the effect of point sources 
leads to a statistical overestimate of the background (despite no 
effect being detected individually). 

We use the 4cr detection threshold in the wavelet analysis 
using the 32" and 64" scales applied to point-source free maps. 
In the following we will report direct flux estimates, made within 
ellipses determined by detection of the signal at 90% confidence 
level. When the signal strength on flux measurement is below 
4cr, simulations sho w that contamination from unresolved point 
sources can be high dBurenin et al.ll2008h . 

3.2. Cluster Red-Sequence 

In order to identify optical counterparts of the extended X- 
ray sources, we h ave made use of the red-sequence technique 
used in .Finoguenov et al.l (1201 Oh . which is its elf a refinement of 
i Jie ph oto-z concentration technique used by IFinoguenov et al.1 
(l200l . 

The method is based on the modeUng of the red-s equence, 
which is performed using the lBruzual & Chariot! (1200 3h popula- 
tion synthesis code, with a passive evolution model of a single 
stellar po pulation ("SSP) and assuming the Chabrier initial mass 
function (IChabrienl2003D and no dust-extinction. The slope of 
the red-sequence is then reproduced by fitting the SSP models 
(formed at Zf = 5) to the Coma cluster red-sequence with a range 
of metallicities. We assume that the slope of the red-sequence is 
entirely due to the mass-metalicity relation, as suggested by both 
observations and theoretical work (IKodama & Arimotd 119971 : 
Stanford etal...l998.) . 

For a given X-ray source, we fit galaxy colours and magni- 
tudes from the K, selected catalogues at the cluster location to 
the red-sequence model described above over the redshift range 
< z < 2. At each redshift step within this range, galaxies within 
an area of radius 0.5Mpc of the centre of the X-ray detection 
and within Az;,/,o/ < 0.2 of the given redshift step are extracted 
and the significance of an overdensity of red galaxies around the 
model red-sequence is estimated. This significance is determined 
with a weighting applied such that galaxies closer to the X-ray 
centre are given a higher weighting. All significant red-sequence 
detections are then recorded. The weighting, w(z), is given by: 



w(z) = 2 



exp 



Ci - Cm(z, nii) 



o-i 



■ '<(2) 






(1) 

where c, and m, are the observed colour and magnitude of the 
/th galaxy, ctj^c is the error on the observed colour, c,„(z, m,) is 
the model red-sequence colour at the observed magnitude, m,, 
whilst ?M*(z) is the characteristic magnitude based on the model 
and cTmag is a smoothing parameter. Finally r, is the distance from 
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Fig. 2. X-ray contour maps of the CFHTLS Dl (left) and D4 (right) fields. The grey-scale shows the smoothed X-ray flux. Detected 
extended X-ray cluster candidates are ringed and numbered, whilst point-source X-ray cluster candidates are marked by red crosses 
and numbered field. Both maps cover an area of 1 ° x 1 °, the entirety of which is covered by the CFHTLS Deep optical u*griz data. 
The sohd fines denote tfie extent of our near-infrared JHKs imaging data. 



the X-ray centre and err is a spatial smoothing factor. The signif- 
icance of the red-sequence measure is estimated based on this 
weighting calculation relative to the median red-sequence signal 
measured at random positions and redshifts within the field. No 
explicit magnitude cut is placed on the galaxy sample for this 
analysis, however the presence of the errors on the colour and 
the magnitude in the denominators means that faint sources with 
high magnitude errors are given low weightings. 

The three different terms in equation \T\ perform weight- 
ing on colour, magnitude and spatial clustering respectively. 
Appropriate colour-magnitude scales are chosen based on the 
redshift range being probed. These are as follows: 

- < z < 0.5: (b - i) colour, i magnitude 

- 0.5 < z < I', (r - z) colour, z magnitude 

- 1 < z < 1.5: (/ - 7) colour, J magnitude 

- 1.5 < z < 2: (z - Ks) colour, Ks magnitude 

As stated, we focus solely on the high redshift sources identi- 
fied by the red-sequence analysis, as it is at z > 1.1 that the deep 
near-infrared data provides significant improvements ov er what 
is po ssible using optical data alone. For example, Ols en et aLl 
(12007 ) perform a cluster search in the CFHTLS fields using the 
CFHTLS deep optical data and are limited to identifying red- 
shifts of counterparts for clusters at z < 1.1 only. 

The success of applying the red-sequence method to iden- 
tifying cluster members at z > 1.1 is reliant on a number of 
factors. The first factor is the sensitivity of the data to the pres- 
ence of the red-sequence. We show the characteristic magni- 
tude as a function of redshift in the top panel of Fig. |3] In the 
range 1 < z < 1.5, this is / band magnitude is used, whilst at 
1.5 < z < 2.0, the Ks band magnitude is used. Additionally, 
we show the median photometric redshift errors as a function of 
redshift for all galaxies within m*, + 1 in both the Dl (solid line) 



and D4 (dashed line) fields (middle panel) and the n(z) distribu- 
tion of this population again for both the Dl (solid line) and D4 
(dashed line) datasets (bottom panel). The increase in the photo- 
metric redshift errors will act to smear the redshift distribution. 
However, for the m = m*„ + 1 population potted, the median pho- 
tometric redshift error remains well below the redshift selection 
window (zrs ± 0.2) we apply prior to the red-sequence analysis. 

Furthermore, the method is reliant on the red-sequence be- 
ing in place in clusters in the survey volume. At the present 
time, the red-sequence has been confirmed to be in place at leas t 
to z = 1.62 based on the observations of iTanaka et al.l (l2010l) . 
However, observations of clusters at z > 1 .5 are somewhat lim- 
ited and the prevalence of the red-sequence and it's dominance 
over star-forming galaxies in clusters at these redshifts is far 
from certain. A fraction of groups and clusters at such redshifts 
may be dominated by star-forming galaxies, in which case the 
success rate of the red-sequence analysis would be reduced. 

We note that the model red sequence shows an off'set with 
respect to the observed red sequence at a given redshift. This is 
due to the systematic uncertainty in the models and also due to 
small photometric zero point errors in the observed data. We do 
not calibrate the off'set with spectroscopically confirmed clusters 
as it is difficult due to the lack of spectroscopic redshifts par- 
ticularly at high-z. The cluster redshifts are estimated using the 
median redshifts of the galaxies identified as cluster members 
via the above method. 

For all of the detections, we perform visual checks of the red- 
sequence identifications and assign a flag based on the quality of 
the identification. These were assigned based on the following 
criteria: 

1 . Good X-ray centre, reliable redshift estimate. 

2. Unable to use X-ray centre, reliable redshift estimate. 

3. Spectroscopic redshift confirmation required. 

4. Potential projection eff'ects. 
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Fig. 3. Top panel: The characteristic magnitude as a function of 
redshift used in equation [T] Note that the J band magnitudes 
are used in the redshift range 1 < z < 1.5 and the K^ band 
magnitudes are used at 1.5 < z < 2.0. Middle panel: Median 
photometric redshift errors of galaxies in the magnitude range 
m*,(z) + 1 as a function of redshift for the Dl (solid line) and 
D4 (dashed line) fields. Bottom panel: The photometric redshift 
distribution of galaxies in the K^ selected catalogue with magni- 
tudes of m - m*J(z) + 1 . 



5. Redshift identification potentially unreliable. 

6. Identification ok, but not a group/cluster. 

7. No identification possible. 

8. Detection out of IR data area. 

Note that as we only present a subset of the cluster identifica- 
tions in this paper, not all the flag values are utilized in this paper, 
however we list the full range of possible flags for consistency 
with the full catalogue, which includes the z < 1 . 1 detections. 



3.3. Cluster Properties & Scaling Relations 

From the cluster apperture flux measurements and redshift esti- 
mates we derive the total cluster fluxes, luminosities and masses. 
The XMM X-ray observations will not be sensitive to the outer 
faint regions of the cluster X-ray emission, particularly so given 
the distant nature of the clusters presented here. Deriving an ac- 
curate estimate of the total flux therefore relies on two possi- 
ble methods, re-observation to gain greater depth or modeling of 
the X-ray profile. In this work we extrapolate the total flux for 
each cluster from the observed flux, Fd, by iteratively fitting a/3- 
model to the observed X-ray emission, from which we estimate 
the correction factor, C^(z, T). Thus the total flux is given by: 



^(< rsQo) = Cpiz, T)Fa 



(2) 



where rjoo is the radius enclosing the matter density 500x the 
critical density, which corresponds to the observed steepening in 



the surface brightness profile s of clusters (IVikhlinin et al.ll2006h . 
iReiprich & Bohringerl (l2002h . amongst others, have shown that 
the bulk of the cluster flux is contained wit hin this radius. The 
j8-model (ICavaliere & Fusco-Femianolll976h is given by: 



S(r)^S(Q) 



1 -t- 



8+1/2 



(3) 



where S (0) is the central flux, r is the cluster radius and we use 
the parameters j6 and rcore given by the scaling relations: 



P = 0.4(A:r/lkeV) 



1/3 



and 



0.63. 



r,„„ = 0.07(^r/lkeV)""V5oo 



(4) 



(5) 



We note that although the y6-model is considered a good de- 
scription of the cluster shape, characterization of the outskirts 
remains uncertain, whilst groups show a larger scatter in their 
X-ray profiles. Further, given the redshift range we are probing 
in this paper, it is important to note that these scaling relations 
have been derived at redshifts of z < 1 . 

In order to estimate the intrinsic cluster luminosity, we ap- 
ply the g -correction following t he ap p roach used by and de- 
scri bed in IReiprich & B o hri ngej (|2002|) . iBohringer et al.1 (|2004|) 
and lFinoguenov et 30 (12007) among others. Therefore, applying 
the A'-correction to account for the temperature and redshift of 
individual clusters, the intrinsic luminosity is given by: 



io.i-2.4keV = ^ndlKiz, T)Cp(z, T)Fd 



(6) 



where di is the luminosity distance to the cluster. 

From the luminosity, we estimate both the mass and tem- 
peratures for a given clu ster based on derived scaling relations 
(iMarkevitch et al.l Il998h . The cluster temperature is therefore 
given by: 

kT = 02 + 6x lo''"^"'^^"'"^'''^^^^-'"^'*''*^^^^' (7) 

where E. is given by: 



E,^ ^JnM(l + z)^ + ^A 



(8) 



The radius measure, r^oo, is then related to the cluster tem- 
perature, kT, by: 



0.63 , 



rsoo = 0.39 IMpc (kT /keWy"' /E, 



(9) 



We solve for the above inter-related parameters iteratively 
for each cluster starting from the observed X-ray signal. 

The above procedure of measuring the cluster luminosity 
has been calibrated against weak lensing mass measurement in 
Leauthaud et al. , (20L0i) within z < 1, yielding 



M200E: = lO'^-^Mo (1.07 + 0.15) 



LyE: 



lO^Z-^ergs.s- 



0.64±0.03 



(10) 



This is consistent with local measurements of the mass- 
luminosity relations, within the adopted coiTection for expected 
evolution with redshift. In this work we perform the mass es- 
timates of the clusters by extrapolating this relation to higher 
redshifts. 
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Table 3. WIRD Cluster Survey z > 1.1 cluster candidates. 



ID 


Field 


Cat no.^ 


R.A. Dec. 

(12000) 


Zrs 


Crs 


Zphot 




Flag 


WIRDXC 10224.3-0408 


Dl 


23 


36.0777 


-4.1403 


1 26+"-" 


9.1 


1.45 





3 


WIRDXCJ0225.0-042I 


Dl 


44 


36.2442 


-4.3550 


1 24+''-2f' 


6.5 


1.39 





3 


WIRDXC J0227.2-0423 


Dl 


48 


36.7988 


-4.3861 


1 24+''-0'* 

^•^^-0.09 


2.9 


1.24 


1.32 


5 


WIRDXC J0225.2-0429 


Dl 


58 


36.2938 


-4.4955 


1 46+'''''* 

^•^"-0.08 


4.0 


1.46 


— 


5 


WIRDXC J0225. 8-0434 


Dl 


64 


36.4403 


-4.5826 


1 1/1+0.04 

^•^^-o.oi 


2.7 


1.14 


1.10 


5 


JKCS 041 (a) 


Dl 


76 


36.6808 


-4.6951 


1 07+0"* 

^•"'-0.06 


5.4 


1.10 


1.13 


1 


JKCS 041 (b) 


Dl 


76 


36.6808 


-4.6951 


1 40+0.28 


6.0 


1.49 


-— 


3 


WIRDXC J2216.3-1729 


D4 


17 


334.0701 


-17.4900 


1 46+"-"'' 


2.1 


1.46 


— 


5 


XMMXCS J2215.9-1738 


D4 


25 


333.9975 


-17.6340 


1 Q7+0.09 


17.9 


1.37 


1.45 


1 


WIRDXC J2216.4-1748 


D4 


32 


334.1084 


-17.8084 


1 41+0.09 
^•^^-0.08 


3.1 


1.40 


1.42 


3 


WIRDXC J2213.9-1750 


D4 


37 


333.4828 


-17.8372 


1 16+"" 


5.4 


1.16 


— 


3 


BLOXJ2215.9-1751.6 


D4 


38 


333.9883 


-17.8574 


1 11+012 


8.5 


1.17 


— 


3 


WIRDXC J2214.2-1757 


D4 


50 


333.5445 


-17.9652 


1 20+"" 


2.5 


1.28 


— 


5 



" Running ID number in full extended X-ray source catalogue 

* Redshifts in italics are tentative based on a single spectroscopic redshift, spectroscopic redshifts of galaxies not selected in the red-sequence 
analysis or the redshift of a QSO. 



4. Cluster Candidates at z > 1.1 

Our XMM-Newton maps of both fields (CFHTLS Dl and D4) 
are shown in Fig. [31 with the z > 1.1 cluster candidates denoted 
by the red ellipses and numbers (the numbers corresponding di- 
rectly to the IDs given in table|4|. The extent of the WIRDS data 
is shown by the solid outlines and corresponds to 5 WIRCam 
pointings in the Dl field and 3 in the D4 field. In total we detect 
62 extended X-ray sources within the Dl field and 40 within the 
D4 field. Typically all sources are based on detections of > 30 
counts. In the Dl, we find that 46 out of the 62 detections appear 
to be associated with galaxy over-densities identified using the 
red-sequence analysis, with 40 being at z < 1.1 and 6 being at 
z > 1.1. In the D4, the number of X-ray sources associated with 
a red-sequence detection is 28 out of the 40, with 22 being at 
z < 1 . 1 and 6 being at z > 1.1. The X-ray detection found to 
not be associated with galaxy over-densities are most likely the 
result of source confusion of close pairs. This is an issue that we 
are actively working to resolve, whilst our requirement of dual 
detection of both an extended X-ray emission and a clustered 
red-sequence should minimize the effect of false detections in 
our final catalogue. As discussed, in this paper we focus on those 
sources with red-sequence identifications at z > 1.1 as it is the 
identification of this population that the deep near-infrared data 
facilitates. 

In table [3] we present the locations of the extended X-ray 
detections and the results of the red-sequence analysis. We list 
the red-sequence redshift estimate, Zrs (column 6) with the sig- 
nificance of the result, cTrs (column 7). The median photometric 
redshift of the galaxies selected via the red-sequence analysis, 
Zphot, is given in column 8. Spectroscopic redshift identifications 
are given in column 9, where redshift in italics represent tenta- 
tive identifications based on either a QSO that may be identified 
with the cluster, non-red-sequence galaxies at the estimated clus- 
ter redshift or just a single red-sequence identified galaxy having 
a spectroscopic redshift. Normal text denotes multiple corrobo- 
rating spectroscopic identifications of red-sequence galaxies. 

For each candidate we present the spatial distribution of red- 
sequence selected galaxies, the magnitude-colour diagram used 



to identify the red-sequence, the red-sequence 'signal' as a func- 
tion of redshift and the redshift histogram within the analysis 
radius. We also provide for each cluster a 2.5' x 2.5' (or larger 
for the more extended dete ctions) gi/T, c olour image constructed 
using the STIFF software (" BertinlUoiC ft. We use the g-band im- 
age to provide the blue channel, the ;-band image for the green 
channel and the ^j-band image to provide the red channel and 
tailor the gamma correction and lower and upper brightness lim- 
its in order to optimize the images for publication. X-ray con- 
tours have been overlayed on the colour images in each case, 
whilst galaxies selected as cluster members via the red-sequence 
method are highlighted (white arrows) and, where available, 
spectroscopic redshifts are denoted (white circles). 



We provide the details of our cluster analysis based on the es- 
timated redshifts and X-ray data in table|4] All candidates listed 
here were detected in the X-ray data and are estimated to be 
at z > 1.1 based on the red-sequence analysis described above 
in conjunction with photometric redshifts. Cluster IDs are given 
(column 1). The best estimated redshift (photometric or spectro- 
scopic) is given in column 2. For each cluster we give r20{) (i.e. 
the radius enclosing a matter density 200x the critical density, 
column 3), the iteratively determined total flux, F{r < rsoo) (col- 
umn 4), and mass, M200 (column 5), the X-ray luminosity, L^ 
(column 6) and the cluster temperature, Tx (column 7). Finally 
the detection reliability flag is given in column 10. Note that for 
the calculations of the cluster properties, we have used a flat cos- 
mology with Qff, = 0.25 and h - 0.72 in order to be consistent 
with die lLeaufliaudetalJ (l2010h . 



All the following candidates were detected using the avail- 
able X-ray data and analyzed using the full 8-band photome- 
try available in the optical from CFHTLS-Deep and the near- 
infrared from WIRDS. We have checked all our candidates 
against known clusters using the NASA Extragalactic Database 
(NED). In cases where a candidate has already been identified, 
with or without a redshift estimate, we compare our results with 
the previous detection. 
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Table 4. WIRD Cluster Survey z > 1.1 cluster candidates. 



ID 


zl 


^■200 


F(< rsoo) 


M200 


L,(0.1-2.4keV) 


Tx 


Flag 






n 


(lO-'^ergs/s/cm^) 


(10'*Mo) 


(10«ergs/s) 


(keV) 




WIRDXC J0224.3-0408 


1.45 


0.0169 


1.92 ±0.31 


0.85 + 0.08 


6.05 + 0.99 


2.18 + 0.15 


3 


WIRDXCJ0225. 0-0421 


1.39 


0.0171 


1.44 + 0.55 


0.73 + 0.16 


4.61 + 2.06 


1.85 + 0.28 


3 


WIRDXC J0227.2-0423 


1.24 


0.0189 


2.43 ± 0.62 


0.91+0.13 


4.96+1.26 


2.10 + 0.22 


5 


WIRDXC J0225. 2-0429 


1.46 


0.0186 


3.02 ± 0.67 


1.03 + 0.13 


7.49+1.67 


2.43 ± 0.22 


5 


WIRDXC J0225.8-0434 


1.14 


0.0213 


4.27 ± 0.84 


1.15 + 0.13 


6.48+1.27 


2.41 + 0.20 


5 


JKCS 041* 


1.13 


0.0259 


< 11.02+ 1.59 


< 1.78 + 0.15 


< 11.68 + 1.68 


3.20 + 0.24 


1 


JKCS 041 


1.49 


0.0224 


< 10.42+ 1.50 


< 1.88 + 0.16 


< 21.42 + 3.08 


3.85 + 0.24 


3 


WIRDXC J2216.3-1729 


1.46 


0.0162 


1.19 + 0.46 


0.68 + 0.15 


3.72+ 1.45 


1.80 + 0.27 


5 


XMMXCSJ2215.9-1738 


1.45 


0.0229 


10.80 + 0.45 


1.91+0.05 


20.50 + 0.85 


3.82 + 0.07 


1 


WIRDXC J2216.4-1748 


1.40 


0.0203 


4.59 + 0.49 


1.25 + 0.08 


9.43+ 1.01 


2.73 + 0.13 


3 


WIRDXC J2213.9-1750 


1.16 


0.0213 


3.65 + 0.74 


1.03 + 0.12 


4.91 + 1.00 


2.17 + 0.18 


3 


BLOXJ2215.9-1751.6 


1.17 


0.0224 


5.13 + 0.39 


1.22 + 0.06 


6.64 + 0.50 


2.47 + 0.08 


3 


WIRDXC J2214.2-1757 


1.28 


0.0207 


4.02 + 0.95 


1.14 + 0.15 


6.83+ 1.62 


2.44 + 0.24 


5 



" Best estimate cluster redshift - photometric redshift unless multiple corroborating spectroscopic redshifts are available. 
* We detect multiple structures with the red-sequence analysis for JKCS 041. However, we are unable to separate the X-ray emissions from the 
different structures and therefore give upper limits to the X-ray properties based on taking the entire X-ray emission for each. 



4.1. CFHTLS D1 Extended Sources 

In the D 1 field we detect 6 z > 1.1 extended X-ray cluster candi- 
dates using the combination of X-ray and NIR data. The X-ray 
data covers the entire Ideg^ of the CFHTLS Deep field, whilst 
the NIR data covers ~ 0.6deg . Our z > 1.1 detections are nec- 
essarily limited to the area covered by the NIR imaging. The 
candidate locations are shown over-layed on the smoothed X-ray 
data in Fig. |2] Coordinates and cluster properties are presented 
in table |4] In the following sub-sections we present each of the 
cluster candidates identified in the CFHTLS Dl field. 



4.1 .1 . WIRDXC J0224.3-0408 (Dl -23) 

In the first high redshift candidate in the Dl field, we find a 
clearly extended X-ray signal (yellow contours in Fig. |4]l with 
a extended flux signal of 6.2cr. The region of the X-ray de- 
tection is somewhat obscured in the optical/NIR data due to a 
bright foreground star We see a strong signal in the red-sequence 
analysis, however the signal is very broad and covers the range 
z ~ 1.1-1.6. Two peaks in the analysis are visible, the first 
at z = 1 .26;^!jAo, with a significance of 9.1 cr and the second at 



z = 1.53^yy4, with a significance of 6.4cr. The broadness of the 
two peaks in the red-sequence analysis suggest that they may be 
caused by the same structure at a single redshift (we note that the 
bright star close to this candidate may be affecting the photom- 
etry). Taking the photometric redshift distribution, we observe a 
large peak at high redshift with z - 1 .45. Given the large uncer- 
tainties on the red-sequence results and the general consensus 
with the photometric redshift peak, we take the median redshift 
of this peak (i.e. z = 1.45) as the cluster redshift. There are no 
apparent red-sequence signals below z ~ 1.1. 

Based on this redshift estimate, we predict a cluster mass of 
M200 - 0.85 X IO'^^Mq and an X-ray luminosity ofLi - 6.1 x 
lO^^ergs/s. 



4.1.2. WIRDXC J0225.0-0421 (D1-44) 

Candidate WIRDXC J0225.0-0421 is shown in Fig. |5] The X- 
ray signal is measured with a significance of ~ 2.6cr. Analyzing 



the red-sequence result, we see a broad peak from z ~ 1 .2 to 
z ~ 1.5 with a significance of ~ 5cr. The primary solution gives 
a redshift of z = 1 .24^qq^. This correlates with a peak in the 
redshift distribution and we find a median photometric redshift 
of the galaxies in this structure of z - 1.39. Looking at the spatial 
distribution of the z ~ 1.39 red-sequence galaxies, we find a 
strongly clustered group within the X-ray contours. A total of 13 
galaxies are classified as red-sequence members at the cluster 
redshift within the observed X-ray contours. Again we see no 
significant signs of structures at redshifts of z = 1.1 using either 
the NIR or optically selected catalogues. 

Taking the z - 1.39 solution, we estimate a cluster mass 
of M200 - 0.73 X IO^'^Mq, an X-ray luminosity of L, - 4.6 x 
lO'^-^erg/s and a cluster radius of raoo - 0.017°. 



4.1.3. WIRDXC J0227.2-0423 (D1-48) 

We find a ~ 4cr signal for this X-ray source, which is clearly ex- 
tended (Fig. |6]l. The red-sequence gives two possible solutions, 
one at z = 1.24^!]!]g with a significance of 2. Oct and the second 



atz 



-0.02 



with a significance of 3.3cr. 



A number of the z ~ 1 .24 red-sequence selected galaxies lie 
within the X-ray signal, whilst several of the brightest lie to the 
North- West of the X-ray detection. One of these has a spectro- 
scopic confirmation of z = 1.325 and lies close to the boundary 
of the X-ray signal. A further galaxy selected as being a clus- 
ter member via the red-sequence method (and also close to the 
edges of the X-ray detection) is found to have spectroscopically 
measured redshift of z = 1.481, somewhat higher than our esti- 
mated cluster redshift. The full spectroscopic sample of galaxies 
within ~ 1 .2' of the X-ray centre is given in table |5j ordered 
by their distance from the centre of the X-ray detection, Ar. We 
also provide our photometric redshift estimate (Zphot) for each 
object where available. The flag value provided gives the confi- 
dence level of the spectroscopic redshift measurement (Zspec) as 
described in section lZSl 

We estimate a cluster mass of M200 = 0.91 x 1O''*M0, an 
X-ray luminosity of L^ = 5.0 x 10'*"'erg/s and a cluster ra- 
dius of r2oo = 0.019° based on our z = 1.32 cluster redshift 
estimate. The high sensitivity of X-rays towards low redshift 
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Fig. 4. Top: Spatial distribution, colour-magnitude plot and red- 
shift distribution/significance from the red sequence analysis. 
The red points in the spatial distribution and colour-magnitude 
plot show the photometrically identified red-sequence galaxies. 
Blue points show galaxies not selected as being part of the red- 
sequence at the estimated cluster redshift and grey points show 
galaxies within the cluster radius at other redshifts. The third 
panel shows the redshift distribution of galaxies within the clus- 
ter radius (divided by the overall field redshift distribution - field 
black histogram), whilst the significance of the red-sequence de- 
tections as a function of redshift is given by the solid red line. 
The dashed horizontal lines give the significance scale. Bottom: 
Colour image combining CFHTLS g, i and WIRDS K, band 
imaging of cluster candidate WIRDXC J0224.3-0408. The yel- 
low contours show the X-ray emission. White arrows denote 
galaxies selected via the red-sequence analysis. The area cov- 
ered by the image is 2.3' x 2.3'. 



groups means that the contribution to the X-ray detection from 
the z=0.23 component will be very important, even in the pres- 
ence of a confirmed high-z cluster. In this case the L^ and M200 
should be treated as upper limits. 



Fig.5. As in Fig.S but for candidate WIRDXC J0225.0-0421. 

Table 5. Spectroscopically confirmed galaxies around WIRDXC 
J0227. 2-0423 close to the predicted cluster redshift. 



R.A. Dec. 

(J2000) 


Ar 

n 


(AB) 


Zphot 


^•spec 


Flag 


36.7946 


-4.3777 


0.0094 


21.09 


1 40+""* 


1.325 


2 


36.8162 


-4.3820 


0.0179 


24.87 


1 17+0.06 
^•^ '-0 05 


1.191 


2 


36.7994 


-4.3842 


0.0020 


21.14 


1 irQ+0.06 
^•■^■-'-0.05 


1.179 


1 


36.7836 


-4.3777 


0.0174 


20.45 


1 iS+oov 


1.298 


1 


36.8152 


-4.3797 


0.0176 


23.42 


1 99+0.06 
^•^ -0.05 


1.293 


3 


36.7828 


-4.3779 


0.0180 


22.87 


1.81!°?« 


1.227 


22 


36.8174 


-43981 


0.0221 


22.99 


— 


1.150 


2 



4.1.4. WIRDXC J0225.2-0429 (D1-58) 

X-ray source 58 the Dl field is a flag 5 cluster candidate with 
an estimated redshift of z = 1.46. The red-sequence analysis re- 
sults and thumbnail of the candidate are shown in Fig. IT] The 
red-sequence detection at z = 1 .46 is relatively weak compared 
to other candidates, but we see few signs of alternative red- 
sequence redshifts based on the analysis. The X-ray signal is 
detected at a level of ~ 4cr and we include this source in the high 
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Fig. 6. As in Fig.S but for candidate WIRDXC J0227.2-0423. Fig. 7. As in Fig.S but for candidate WIRDXC J0225.2-0429. 



redshift catalogue to help in indicating the upper limits on the 
numbers of high redshift clusters in this field. 

Iteratively solving for the cluster properties using the scaling 
relations and the z - 1.46 redshift solution, we find a cluster 
mass of M2()o = 1.03 + 0.13 x 10'"* Mq, an X-ray luminosity of 
L,. = 7.49 + 2.43 x 10"*^^ ergs/s and a radius of raoo = 0.0186°. 

4.1.5. WIRDXC J0225.8-0434(D1 -64) 

The thumbnail for candidate WIRDXC J0225.8-0434 is shown 
in Fig. [8] The extended X-ray signal, detected at a signal of ~ 5cr, 
shows a relatively compact structure. The red-sequence analysis 
estimates a cluster redshift of z = 1.13+^^5 with a confidence of 
~ 2.60-. 

We find three Zphot ~ 1-13 red-sequence selected objects 
within the X-ray contours, with a number of other sources at 
comparable photometric redshifts close-by. A number of spec- 
troscopically observed objects are available around this object 
(see Fig.|6ll, two of which lie close to the estimated cluster red- 
shift: the first at Zspec = 1.169 and the second at Zspec - 1-123. 
We note that neither is selected as being a cluster member via 
the red-sequence analysis. 



Table 6. Spectroscopically confirmed objects around WIRDXC 
J0225. 8-0434 within Az = 0.1 of the estimated cluster redshift. 



R.A. Dec. 

(J2000) 


Ar 




(AB) 


Zphot 


Zspec 


Flag 


36.4405 


-4.5790 


0.0036 


23.81 


1 06+""' 

^•""-0.06 


1.169 


3 


36.4332 


-4.5798 


0.0076 


22.61 


1 07+""3 

^•"'-0 03 


1.123 


3 


36.4377 


-4.5814 


0.0029 


20.69 


1 nc+0.04 
i-UJ_0.04 


1.095 


2 


36.4383 


-4.5763 


0.0066 


21.73 


1 19+0.03 


1.176 


3 


36.4432 


-4.5748 


0.0083 


21.05 


1 9Q+0.03 
^•^■^-0.03 


1.177 


3 


36.4228 


-4.5884 


0.0185 


22.65 


1 19+0.04 
^•^^-0.03 


1.170 


2 



Taking the estimated cluster redshift of z ~ 1.13, we estimate 
a cluster mass of M200 = 1.15 x lO^'^M©, an X-ray luminosity of 
Lv - 6.5 X lO'^-'erg/s and a cluster radius of r2oo - 0.021°. 

4.1.6. JKCS041 (D1-76) 

This candidate was reported by JAndreon et al.l (l2009l) . It is 
shown in Fig. [10] and appears to be a complex combination of 
different structures along the line of sight. The X-ray detection is 
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Fig. 8. As in Fig.gl but for candidate WIRDXC J0225.8-0434. 
Spectroscopic redshifts are shown in green in the lower panel. 



clearly extended with a detection of 1 0.4+ 1 .5 x 1 0" '^ergs/s /cm^, 
making it the second brightest object in ourz > 1.1 cluster cat- 
alogue. Within the extended X-ray emission we also detected a 
point-source signal in the XMM data, the location of which is 
given by the blue x in Fig. [TO] This point-source was identified 
and filtered out by the point-source removal processing and was 
measured to have a flux of 5.6 x 10"'^ergs/s/cm^. 

From the public V VDS Deep Survey, a large number of spec- 
troscopic redshifts are available within the field. Analyzing the 
distribution of the spectroscopic identifications within 1.2' of the 
X-ray profile centre (black histogram in Fig. |9]l, three redshift 
peaks are evident at redshifts of z = 0.80, z - 0.96 and z - 1.13 
(the redshift distribution of the complete VVDS deep sample is 
given by the filled grey histogram in Fig.|9]l- We plot the galaxies 
identified at these redshifts in Fig.[TO] Identified z ~ 0.8 galaxies 
are given by green circles, z ~ 0.96 galaxies by orange circles 
and z~ 1.12 galaxies by red circles. 

Combining this data with the photometric redshift data and 
the red-sequence analysis, we find firstly that the z - 0.80 galax- 
ies appear to be offset from the X-ray emission, whilst we see 
no significant signal in the red-sequence analysis. From the pho- 
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Fig. 9. Spectroscopic redshift distribution within 1 .2' of the cen- 
tre of the X-ray emission for candidate JKCS 041 (solid black 
line histogram). The redshift distribution of the entire sample is 
shown by the filled grey histogram. 



tometric redshift analysis, we find that the photometric data for 
these galaxies are (where an accurate photometric redshift has 
been calculated) best fit by starburst templates, corroborating the 
lack of any red-sequence in these objects. Given the spatial dis- 
tribution of these galaxies and their blue colours, it is likely that 
these are not associated with the bulk of the extended X-ray sig- 
nal. 

Taking the z - 0.96 spectroscopic redshift peak, we find a 
small signal in the red-sequence analysis, whilst the positions of 
the spectroscopically confirmed z = 0.96 galaxies are well cor- 
related with the X-ray contours. Again however, for the galax- 
ies for which the photometric redshifts corroborate the spectro- 
scopic redshifts, the photometric data are best fit by starburst 
templates. Again the blue colours of these galaxies suggest that 
they may not be associated with the bulk of the X-ray produc- 
ing gas, although they remain a strong candidate for some sort 
of structure, potentially a small group or filamentary structure 
along the line of sight. 

Moving to the third solution, we find a strong signal from the 
red-sequence analysis at z = 1 .07 + 0.06, which corroborates the 
z - 1.13 estimate from the spectroscopic redshift distribution. 
Taking the photometric fitting, the photometry for two of the 
confirmed z - 1.13 galaxies is best fit by elliptical templates. Out 
of the three spectroscopically identified solutions, we take the 
z = 1.13 solution as the most plausible structure to be associated 
with the bulk of the X-ray emission. 

From the red-sequence analysis however, we find the 
strongest red-sequence solution at a redshift of z = 1 ■49^gg^. 
Selected cluster members at this redshift are highlighted by 
white arrows in Fig. [10] and are observed to align well with the 
X-ray contours. Three spectroscopic redshifts at z ~ 1.5 are 
available from the VVDS data, however these are not selected 
as part of the red-sequence analysis and are attributed with spec- 
troscopic flags of only 1 and 2. Despite the lack of spectroscopic 
confirmation, this apparent z - 1 .49 structure is a strong candi- 
date as contributing source to the X-ray detection. 
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Table 7. Spectroscopically confirmed objects around JKCS 041 
within Az = 0.1 of the estimated cluster redshift. 



Cluster ID : Dl-76 



R.A. Dec. 

(J2000) 


n 


(AB) 


Zphot 


Zspec 


Flag 


36.6917 


-4.6949 


0.0108 


20.4 


76+'"" 

"•'"-0.03 


0.797 


4 


36.6920 


-4.6920 


0.0116 


23.1 


n 78+0.04 


0.794 


9 


36.6861 


-4.7099 


0.0158 


19.6 




0.795 


4 


36.6720 


-4.6811 


0.0166 


23.0 


75+004 

"■'-'-0.04 


0.797 


2 


36.6973 


-4.6907 


0.0170 


21.8 


75+003 

"•'-^-0.03 


0.798 


2 


36.6831 


-4.6980 


0.0037 


— 


2.10!«;^t 


0.959 


3 


36.6851 


-4.6883 


0.0080 


22.1 


0.87!«- 


0.960 


3 


36.6848 


-4.6870 


0.0090 


20.4 


QQ+0.04 
"•"-0.03 


0.962 


2 


36.6715 


-4.6930 


0.0095 


21.2 


1 2n+""3 


0.963 


2 


36.6713 


-4.6832 


0.0153 


23.0 


n Q4+0.04 
"•^^-0.06 


0.962 


2 


36.6779 


-4.6763 


0.0190 


20.8 


0.82!««^ 


0.965 


2 


36.6822 


-4.6978 


0.0030 


24.8 


0.95t]? 


1.127 


22 


36.6751 


-4.7030 


0.0098 


20.3 


1 04+0.03 


1.128 


2 


36.6687 


-4.6911 


0.0128 


21.6 


1 n8+0.03 


1.130 


3 


36.6674 


-4.7059 


0.0173 


20.0 


1 12+0'"' 

^•^^-0.04 


1.125 


3 


36.6648 


-4.7045 


0.0185 


19.6 


96+'"'3 

"•^°-0.03 


1.125 


3 


36.6951 


-4.6966 


0.0143 


24.0 


T 78+0.04 


1.490 


2 


36.6891 


-4.6798 


0.0174 


22.4 


1 56+"''* 


1.517 


2 


36.6996 


-4.6890 


0.0198 


22.2 


1 77+0.13 

^•"-0.12 


1.537 


1 



As stated, this X-ray source is at the same position as that 
reported by Andreon et al. (2009) as a z = 1.9 cluster (the cen- 
tre coordinat es of the two detec tions are within 0.2' of each 
other). The Andreon et al.l (l2009h X-ray detection was made us- 
ing observations from the Chandra X-ray telescope, making this 
XMM based observation an independent detection of this clus- 
ter Andreon et al. (2009) performed a red-sequence analysis us- 
ing UKIDSS DR3 deep J and K NIR data. The latest depths for 
these bands re ported for the UKIDSS Seep Extragalactic Survey 
dWarren et al.ll2007h are /(Vega) = 22.2 and 7i:(Vega) = 20.9 
given 5cr detected point sources with 2" apertures (~ 23.1 and 
~ 22.8 respectively in the AB system). This compares to depths 
(in the AB system) from the WIRDS data of 7 = 24.3, // = 24.1 
and ^.5 = 24.1, again using 5cr detected point sources with 2" 
aperture. 

Using the UKIDSS JK data. lAndreon et al.1 (l2009l) reported 
a cluster redshift of z = 1.9 with a confidence of 6.5cr. In con- 
junction with this detection, they also comment that at least one 
more cluster is detected along the line of sight. Based on the 
red-sequence analysis with the CFHTLS/WIRDS data, we find 
no significant detection of a red-sequence at z = 1.9. 

We make estimates for cluster properties taking the z = 1.13 
andz = 1.49 candidate structures. As discussed, the z = 1.13 so- 
lution appears to be confirmed by spectroscopic observations, 
whilst the z = 1.49 solution remains tentative and based on 
photometric/red-sequence observations. Taking these two solu- 
tions, and using the total extended flux measurement in each 
case, we place upper limits on the cluster mass for a cluster at 
z = 1.13 of M200 < 1-78 + 0.15 X IO'^^Mq and for a cluster at 
z= 1.49ofM2oo < 1.88+O.16x1O14M0. We note that based on 
this estimated upper limit for the mass of the possible z - 1 .49 
structure that the spectroscopic redshifts of the three objects span 
too large a range to all be part of the same bound structure. 
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Fig. 10. As in Fig.|4l but for candidate JKCS 041 and with an im- 
age size of 2.8 ' x2.8' . The coordinates of the cluster as given by 
And reon et al.l (l2009l) are shown by the blue cross. The location 
of the filtered out point-source is given by the blue x. 



4.2. CFHTLS D4 Extended X-ray Sources 

We have a total of six z > 1.1 extended X-ray cluster candidates 
identified in the CFHTLS D4 field based on the initial X-ray de- 
tections in combination with the red-sequence/photometric red- 
shift analysis. The distribution of these candidates are shown 
(circled in red) in Fig.|2] over-layed on the filtered X-ray map. 



4.2.1 . WIRDXC J221 6.3-1 729 (D4-1 7) 

The X-ray detection for candidate WIRDXC J2216.3-1729 is de- 
tected at a level of ~ Icr. This is a relatively weak source, with 
a tentative red-sequence detection at z = 1.46. This detection (at 
a level of ~ 2.5cr) correlates with a peak in the photometric red- 
shift distribution at the same redshift. Given the low detection 
signals involved, we attribute this candidate flag 5 and include it 
in the sample as part of placing upper limits on the numbers of 
clusters in the fields. 
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Fig. 11. As in Fig.S but for candidate WIRDXC J2216.3-1729. 



From the measured X-ray flux and the redshift estimate pre- 
sented here, we find a cluster mass of M2()o = 0.68 +0.15 x 
IO'^^Mq, an X-ray luminosity of L, = 3.72+ 1.45 x 10'*^ ergs/s 
and a radius of raoo = 0.0162°. 



4.2.2. XMMXCS J221 5.9+1 738 (D4-25) 

This detection is the same as the spe ctroscopically confirmed 
z = 1.4 cluster of IStanford et al.l ( 1200 6). From our imaging, we 
see a strong clustering of red-sequence galaxies both spatially 
(Fig. [12] lower panel) and in magnitude-colour space (Fig. fT2l 
upper panels). The red-sequence analysis and photometric red- 
shift distribution both indicate a cluster redshift o fz = 1.37. This 
is rea ssuringly close to the redshift measured by IStanford et al.l 
( 1200 6) of z = 1 .45 and confirms the reliability of our redshift 
estimation techniques. In Fig.[T2] we show those galaxies iden- 
tified as cluster members by our red-sequence analysis. Four of 
th ese have been spectro scopically confirmed as cluster members 
by Stanford et al. (2006) and are circled and labeled in the figure. 
[Stanford et al. (2006) present spectra for a further two cluster 
members that are not selected as cluster members as part of our 
analysis, whilst from the AAOmega spectroscopic data we ob- 



Fig. 12. As in Fig.g] but for candidate XMMXCS J2215.9-1738 
and with an image size of 3.3' x 3.3'. The white circles and 
numbers in the lower panel indicate objects with spectroscopic 
redshifts. The blue cross sho ws the coordinates of XMMXCS 
J2215.9-1738 as published byHtanfordeFiDdlOOi). 



serve a QSO at z = 1 .462 lying in what appear to be the cluster 
outskirts. We give all the available spectroscopic redshifts from 
[Stanford et al. (2006) and AAOmega in table[H] 

Based on the spectroscopically measured cluster redshift 
of z = 1.45, we estimate a cluster mass of M200 = 1.9 + 
0.05 X 10'"* M0. Taking the velocity dispersion of cr,, - 580 ± 
140 kms ', measured bv ' Hilton et al.l (l2007l) . and combining 
this with the virial radius of R^, - 1 .05 Mpc estimated by the 
same authors, the total virial mass is estimated to be Mtoi = 
1.5 + 0.6 X lO^'* Mq, comparable to our estimate from the X- 
ray observation. The X-ray luminosity is estimated to be L^ = 
20.50+0.85 xlO"*^ ergs/s, with a cluster radius of r2()o = 0.0229°. 



4.2.3. WIRDXC J221 6.4-1 748 (D4-32) 

WIRDXC J2216.4-1748 is a clear- extended X-ray detection, 
with several red-sequence galaxies at z = 1 .4 found within the 
X-ray detection (Fig.flJl). The significance of the z = 1.4 red- 



Table 8. Spectroscopically confirmed objects around XMMXCS 
J2215.9-1738 within Az = 0.1 of the estimated cluster redshift. 
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Table 9. Spectroscopically confirmed objects around WIRDXC 
J2216.4-1748 within Az - 0.1 of the estimated cluster redshift. 
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sequence detection is 2.7cr, whilst the X-ray detection is mea- 
sured with a signal of ~ 9cr. Although the red-sequence analysis 
does not appear conclusive, we note that the detection coincides 
with a strong peak in the photometric redshift distribution and 
that no significant signal in the red-sequence analysis is seen at 
any other redshifts. 

We note that, although we have no spectroscopic confirma- 
tion of selected cluster members, we observe a spectroscopically 
confirmed z - 1 .42 quasar (taken from the AAOmega data and 
given in table |9]l within the extended X-ray detection, which we 
speculate may be part of this structure. 

For this candidate, we find a cluster mass of M200 - 1 -25 + 
0.08 X 10'^ Mq, an X-ray luminosity of L,- = 9.43 + 1.01 x 
10'*-' ergs/s and a radius of raoo - 0.0203°. 

4.2.4. WIRDXC J221 3.9-1 750 (D4-37) 

The red-sequence analysis and giK^ image with selected cluster 
members and spectroscopic identifications for cluster candidate 
WIRDXC J2213.9-1750 are shown in Fig.[Il The X-ray detec- 
tion for this source is measured at a signal of ~ 5cr and is clearly 
extended. 

The detection lies close to the edge of our NIR data and 
thus have larger photometric errors than in the main regions of 
our data. In addition the central region is obscured by a bright 
star, hampering the red-sequence analysis. Despite these prob- 
lems, we detect a strong peak in the red-sequence analysis at 
z - 1.16 with a significance of 4.6cr. The galaxies selected as 
cluster members are well clustered within the extended X-ray 
emission. 

We have only a single spectroscopic redshift from the 
AAOmega data in the region of this candidate, however it lies 
at a lower redshift and is not selected as a cluster member by our 
analysis. Based on a cluster redshift of z = 1.40, we estimate a 
cluster mass of M200 = 1.0 x lO'^^M©, an X-ray luminosity of 
L, = 4.9 x lO'^^^erg/s and a cluster radius of r2oo = 0.021° 

4.2.5. BLOX J2215.9-1751.6 (D4-38) 

This c andidate is also detected in the work of iDietrich et alJ 
( l2007l) and is listed as BLOX J2215.9-1751.6. They find a re- 
liable X-ray detection with a flux measurement of 5.7 + 0.6 x 
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Fig. 13. As in Fig.lH but for candidate WIRDXC J2216.4-1748 
and with an image size of 2.5' x 2.5'.. The white circles and 
numbers in the lower panel indicate objects with spectroscopic 
redshifts. 



10 '^ ergs/s/cm^ and an angular diameter of the major axis of 
16.8". This compares to our detection ofF{< rsoo) = 5.1 +0.59x 
10"'^ ergs/s/cm^ and r2oo = 0.0224° {^ 80"). They are unable 
to measure a redshift for this source however as their analysis 
is limited to BVRI imaging with no NIR data. The candidat e 
may also be present in the cluster survey of lOIsen et al.l (l2007h . 
who present a candidate (CFHTLS-CL J221556-175023) at a 
distance of 1.05' from o ur detected c luster core with a redshift 
estimate of z = 1.1. The lOlsen et alJ (2007 ) analysis is limited 
to CFHTLS u*griz data and attribute the detection their lowest 
confidence grade (D). The locations of the X-ray detections for 
BLOX J2215.9-1751.6 and CFHTLS-CL J221556-175023 are 
given by the blue crosses in the giK,, image in Fig. [15] BLOX 
J2215.9-1751.6 aligns with the centre of our X-ray contours, 
whilst CFHTLS-CL J22 1556- 175023 is offset to the North. 

Using the optical plus NIR data, our red-sequence analysis 
(Fig. \T5[ produces a significant peak in the probability distribu- 
tion at z = 1.11 (with a significance of 8.5cr). This is reinforced 
by a strong peak in the photometric redshift distribution with a 
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Fig. 14. As in Fig.g] but for candidate WIRDXC J2213.9-1750 
and with an image size of 2.5' x 2.5'. 



median redshift of z = 1.17. Galaxies selected as part of the red- 
sequence at z = 1.17 are highlighted with the white arrows in 
Fig. [15] A number of the brightest of these galaxies are clearly 
clustered within the core of the X-ray extended signal and appear 
to follow the shape of the X-ray profile. 

A small number of spectroscopic redshifts are available in 
the region around this cluster candidate, however none are avail- 
able for any of the galaxies selected as cluster members via the 
red-sequence analysis. 

Estimating the clustering properties from the given redshift, 
we find a cluster mass of M200 - \.T1± 0.06 x 10'"* Mq, an 
X-ray luminosity of L^ = 6.64 + 0.50 x lO"*"* ergs/s and a radius 
of rzoo = 0.0224°. 

4.2.6. WIRDXC J221 4.2-1 757 (D4-50) 

The thumbnail for candidate WIRDXC J2214.2-1757 is shown 
in Fig. [16] The extended X-ray signal, detected at a signal of 
~ 4cr, shows a relatively compact structure. The red-sequence 
analysis estimates a cluster redshift of z = 1 .28 with a confidence 



Fig. 15. As in Fig. |2 but for candidate BLOX J2215.9-1751.6 
and with an image size of 3.3' x 3.3'. The blue crosses mark the 
ce ntral coordinates of th e detections of BLOX J2215. 9-1751. 6 
bv iDietrich et~an (12007) and of CFHTLS-CL J221 556-175023 

ofl'set by 1.05' northwards of our detection) by lOlsen et al] 

I2007h . 



Table 10. Spectroscopically confirmed objects around WIRDXC 
J2214. 2-1757 within Az = 0.1 of the estimated cluster redshift. 



R.A. Dec. 

(J2000) 



Ar 

(°) 



(AB) 



^pho 



Class 



333.5403 -17.9484 0.0173 19.38 



1.196 QSO 



of 3cr. Visually the red-sequence selected galaxies do not appear 
to be well spatially clustered. 

We have no spectroscopic data for any of the galaxies in the 
X-ray emission region. However we note another QSO from the 
AAOmega spectroscopic data close to the X-ray emission and at 
a redshift of z = 1.196. This is given in table [TOl 

Based on a redshift of 1.28, we estimate a cluster mass of 
M200 - 1-14 + 0.15 X 10'"* M0, an X-ray luminosity of L ^ = 
6.83 + 1.62 X lO''^ ergs/s and a radius of raoo = 0.0207°. 
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Fig. 16. As in Fig.i] but for candidate WIRDXC J2214.2-1757. 
The white circles and numbers in the lower panel indicate ob- 
jects with spectroscopic redshifts. 



4.3. X-ray Unresolved Candidates? 

Thus far all the detected candidates have been part of the ex- 
tended source catalogue. As discussed, this excludes any sources 
that are close in size to the PSF of the image within the con- 
straints we have set. Especially as we are concerned with high- 
redshift systems, this risks the possibility of missing particularly 
compact/faint clusters. We have therefore reviewed the point- 
source catalogue for the two fields, running the red-sequence al- 
gorithm with the weighting tuned to select more compact spatial 
clustering by setting (Tr = 300kpc. 

From this analysis, we find three candidates for high redshift 
compact clusters that are filtered out by the algorithm used to de- 
tect extended X-ray emission. The thumbnails and red-sequence 
results for these are shown in fi gures [TTl [TS] and [T9] 

The first candidate is found in the Dl field and based on the 
red-sequence analysis for this source, we estimate a redshift of 
z = 1.35. This red-sequence detection is found with a strong 
signal, but the red-sequence members are distributed outside of 
the X-ray point source signal. We see no obvious AGN candi- 
dates at the position of the X-ray signal and assign the cluster 



Fig. 17. As in Fig.gl but for candidate WIRDS-CS-Dl-p282. 



candidate a flag of 5. We note that we find a small number of 
objects with spectroscopic redshifts close to the X-ray emission, 
however none of the red-sequence selected objects have been 
observed and the objects that have been observed are not at con- 
sistent redshifts. 

The second candidate found in the point-source list is in the 
D4 field and is shown in Fig. [18] Note that the elongated X-ray 
signal is due to the non-symmetric nature of the PSF across the 
field of view of the observations, however despite this the detec- 
tion is still classed as a point-source by the X-ray detection algo- 
rithm. A relatively strong red-sequence signal is seen, dominated 
by the clustering of red objects to the top-left of the X-ray emis- 
sion. These are found to be at a redshift of z = 1 .26^^ 22 based on 
the red-sequence analysis. We note the presence of a relatively 
bright foreground galaxy which may hold an AGN and therefore 
may be the source of the X-ray emission. Again we attribute this 
candidate a flag of 5. 

The final point-source detection with signs of a red-sequence 
detection is shown in Fig. [19] Again a strong red-sequence signal 
is observed, in this case at a redshift of z = LO^^gg^. A number 
of the selected galaxies are found with the PSF profile of the 
X-ray detection and appear well clustered. A second grouping 
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Table 11. WIRD Cluster Survey z > 1.1 cluster candidates. 



ID 



R.A. Dec. 

(J2000) 



^■200 

n 



F« rioot 



M200 L,(0.1-2.4keV) Tx 

(IO'-^Mq) (lO'^^ergs/s) (keV) 



" Cluster redshift estimated from red-sequence analysis or spectroscopic redshift if available. 
* Extrapolated cluster flux (10"'''ergs/s/cm^). 



Flag 



WIRDCS-l-p282 36.644929 -4.1936979 1.35 0.0162 0.93 ± 0.63 0.60 ±0.22 2.60 ±1.76 1.59 ±0.39 5 

WIRDCS-4-p265 333.9341 -17.699511 1.26 0.0165 0.87 ± 0.20 0.56 ± 0.08 2.05 ± 0.48 1.47 ±0.13 5 

WIRDCS-4-p434 333.3954 -17.872525 1.13 0.0187 1.49 ±0.92 0.67 ± 0.22 2.27 ±1.40 1.57 ±0.35 5 
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Fig. 18. As in Fig.H but for candidate WIRDS-CS-D4-265. 



Fig. 19. As in Fig.H but for candidate WIRDS-CS-D4-434. 



of galaxies with the same red-sequence redshift is seen to the 
upper left of the X-ray emission. Again, a relatively bright fore- 
ground galaxy is observed within the PSF providing a candidate 
for AGN emission as the source of the detection. This candidate 
is also attributed with a flag of 5. 

As with the primary cluster candidates, we estimate the to- 
tal flux/luminosity /mass properties of the compact cluster can- 
didates based on the derived redshifts. The results are given in 
table [TT] We note that these are all low-flux detections at the 



limit of our detection thresholds. As such we estimate relatively 
low-masses for these objects (~ 0.5 - 0.6 x lO'^^M©). 

In summary of the check for unresolved clusters, no obvious 
candidate has been found and the 3 candidates considered are 
likely AGNs and not unresolved clusters. Thus, we can exclude 
a large influence of spatial resolution of XMM-Newton on the 
reported cluster statistics. 
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Fig. 20. Cluster luminosity probed as a function of redshift in the 
Dl (top panel) and D4 (lower panel) fields. The flag=5 (open 
circles) and flag<3 (filled circles) extended cluster candidates 
are shown. Open stars show the point-source confusion cluster 
candidates. The dashed line shows the X-ray luminosity limit 
based on the minimum extended X-ray flux detection in the data, 
whilst the solid line shows the X-ray luminosity limit based on 
the median X-ray flux limit. 



5. The Cluster Sample 

5.1. Cosmology from Cluster Number Counts 

In Fig.|20l we show the derived X-ray luminosities of our clus- 
ter candidates in the CFHTLS Dl and D4 fields as a function of 
redshift. The high-redshift flag<3 cluster candidates are shown 
by filled circles, the flag=5 cluster candidates are shown by the 
open circles and the potentially unresolved detections are shown 
by open stars. For both fields, the dashed line shows the luminos- 
ity limit estimated from the minimum detected extended X-ray 
flux in each field, whilst the solid line gives the luminosity limit 
calculated from the median X-ray flux limit. The point-source 
detections are consistent with being faint sources in which ex- 
tended signal from the cluster outskirts is below the detection 
threshold of the XMM data. 

We present the number counts of clusters as a function of 
redshift in the two fields in Fig. [21] The filled circles show the 
number counts based on only flag<3 candidates and the open cir- 
cles show the cluster number counts for the full flag<5 candidate 
sample. As they include a number of tentative detections, the 
flag<5 counts give the upper limits to the number of detectable 
clusters based on the data used here. In both cases the plotted er- 
rors are statistical. We note that we do not detect any clusters at 
z > 1 .5 and so the final data-point represents an estimated upper 
limit to the cluster number counts based on the areal coverage of 
our survey. 

Model n umber counts w ere predicted using the cosmolog- 
ical code of iPeacocM (12007 ) with th e WMAP 7 -year cos mol- 
ogy (i.e. Q„,/j2 = 0.1334, erg = 0.801. iKomatsu et al.ll2010l) and 
applying the areal coverage and flux limits of our survey. This 
is shown by the solid line in Fig. [21] We find that the model 
over-predicts the number counts compared to the flag<3 candi- 
dates by ~ 2cr, whilst matching well with the full z > 1 . 1 upper 



Fig. 21. Cluster number counts as a function of redshift. The 
filled black circles show the candidate cluster counts from the 
flag<3 objects in the Dl and D4 fields combined. Open cir- 
cles show the counts based on the flag<5 objects. The solid line 
shows the model predictions based on the WMAP 7-year cosmo- 
logical parameters. The effect of the Icr uncertainty in the scal- 
ing relations is illustrated by the dashed line, whilst the dash-dot 
line shows the effect of the Icr uncertainty in the scaling rela- 
tions, combined with a WMAP 7-yr model taking the 2cr uncer- 
tainty limits on Q„, and cr^. 



limit flag<5 sample. There are two key elements that efifect the 
number count measurements and models: the cluster scaling re- 
lations and the cosmological parameters used (in particular Q,„ 
and cTg,). Looking at the effect of changes in both of these, the 
dashed line shows the result of adjusting the scaling parameters 
by Icr, whilst the dash-dot line shows the effect of adjusting the 
scaling parameters again by Icr combined with a change in the 
WMAP 7-year Q.^ and os values by 2cr. The second of these, 
which takes values of Q„,/i^ = 0.1323 and cr^ - 0.741, lowers 
the model number counts by a factor of » 2.5 at redshifts of 
z ~ 1-2.0. This illustrates the sensitivity of the cluster number 
counts to the cosmological parameters and is more consistent 
with our flag<3 candidate counts. The number counts based on 
the flag<3 favour a lower os, which we note is in agreement with 
the results of Finoguenov et al. (2010). Based on their sample of 
< z < 2 clusters, they publish cluster number counts that prefer 
a 5% reduction in the value of cr^ from the WMAP 5-year value. 
Significant uncertainty remains in the number count analy- 
sis. In particular the scaling relations present problems in any 
potential application of cluster number counts to placing con- 
straints on cosmological parameters. The scatter in the scaling 
relations, particularly for galaxy groups, remains significant, the 
effect of which is illustrated by the dashed line in Fig. |2T| In 
addition, we note that the redshift range over which scaling re- 
lations have been cahb rated is in the main limited to z < 1 
(iLeauthaud et al.ll2010t) . Of course, we also note that, although 
we present a relatively large area for such a high-redshift sur- 
vey, statistical errors due to low numbers of clusters remain high. 
Finally, uncertainty in a number of the cluster detections require 
spectroscopic confirmation of cluster members. 
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Returning to the lack of cluster detections at z > 1.5, we 
note that given the WMAP 7-year model we would expect 
«i 3 - 4 cluster detections within the survey area at such red- 
shifts. Based on the red-sequence model, we estimate a charac- 
teristic Ks band magnitude of cluster members of m*^ x 21.5 at 
1 .5 < z < 2, which is well within our magnitude completeness 
limit. Increasing photometric errors will give poorer photomet- 
ric redshifts (see Fig. [3]l, which will weaken the signal from a 
given cluster (note that the photometric redshifts are used to pre- 
select potential cluster members). For galaxies at this character- 
istic magnitude and in this redshift range, we find a median error 
on the photometric redshift measurements of 0.05 < cr, < 0.1, 
compared to « 0.03 - 0.05 at z < 1.5. Looking at the photome- 
try for this same population, we find increases in the magnitude 
errors of ~ 0.05 -O.OSmag. In order to gauge the effect of this in- 
crease in the magnitude errors, we re-run the red-sequence anal- 
ysis on the same catalogues but with a 0. 1 mag error added in 
quadrature to all 8 bands. Taking all the X-ray extended sources 
(including those originally classed as z < 1), we find that 70% 
of detections retain the same redshift solution given this increase 
in the noise of the photometry. Out of the 30% with differing so- 
lutions given the increased photometric errors, 72% have a red- 
sequence significance of cTrs < 3. Based on this analysis then, 
we may expect that the red-sequence completeness at z > 1 .5 is 
reduced by no more than 30%. 

Alternatively, the lack of detections may be the result of 
greater levels of star formation in clusters at z > 1.5. Increased 
star-formation in cluster member may cause the red-sequence to 
be more difficult to detect, if higher fractions of bluer galaxies 
are present. This effect is very difficult to quantify at this point 
as it relies on having samples of unbiased clusters and groups 
at these redshifts in order to understand the properties of clus- 
ter members. This will be further complicated by the formation 
epoch of the red sequence varying from cluster to cluster as clus- 
ters will not all form at the same time. Ultimately, the redshifts 
up to which we can detect clusters via the red sequence method 
is something that must be learned from real data and as such 
the highest redshift at which the red-sequence analysis can be 
claim ed to be effective is z ~ 1.6 at this point (Tana kaet al.l 
120101) . 

Finally, given the low numbers of clusters predicted to be 
detectable by the X-ray data in this redshift range, it is feasible 
that the lack of detection of any clusters at such redshifts may 
be the result of statistical cosmic variance (i.e. that there are no 
clusters above our detection thresholds at these redshifts in the 
volumes sampled). 

5.2. X-ray Luminosity Function 

It is intructive to see in which way the prediction of high num- 
bers of clusters in the WMAP7 cosmology disagrees with ob- 
served cluster characteristics, such as total mass or luminosity. 
In this vein, we calculate the X-ray luminosity function in the 
redshift range 1.1 < z < 2 based on our cluster candidate sample 
in the two analyzed fields. The cluster luminosity function, 0, is 
given by: 
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where A^^ is the number of clusters in a given luminosity bin of 
width AL and ymax(ix,i) is the total comoving volume in which 
a cluster of luminosity L^,, could have been detected above the 
flux limits of the survey. This is given by: 



where Zmin and Zmax are the redshift limits. The lower limit is 
taken as z - 1.1, whilst we set an upper limit of z = 2.0, assum- 
ing for the purpose of the luminosity function that the lack of 
detections at 1.5 < z < 2.0 is due to statistical cosmic variance. 
Q(/x(Lx, z)) is the sky area as a function of the X-ray flux probed 
by the X-ray data and dV{z) ldz is the differential comoving vol- 
ume element per steradian dAvni & Bahcall|[l980l: iMuUis et al.l 
2004. 

Fig. |22] shows our calculated X-ray luminosity function in 
the redshift range 1.1 < z < 2. The luminosity function for the 
flag<3 sample is again shown by the filled circles, whilst the lu- 
minosity function derived from the flag<5 is shown by the open 
circles. Again, errors on the bin numbers are statistical (although 
the horizontal bars simply show the bin extent). The solid line 
gives the model luminosity function based on the WMAP 7-year 
cosmological parameters, whilst the dashed line shows the ef- 
fect of reducing the WMAP 7-year Q^ and erg parameters by 
2cr. We see that the flag<3 luminosity function agrees well with 
the WMAP 7-year model at luminosities of Lx > 4 x 10'*-'ergs/s. 
Taking the single luminosity bin at Lx < 5 x lO^'^ergs/s, we find 
a deficit of clusters compared to the model luminosity function, 
with the cluster data point falling a factor of 3 below the model 
prediction (equivalent to « 2cr based on the error estimate). With 
the flag<5 luminosity function however, we find good consis- 
tency with the model over the full range of luminosities that 
our cluster candidates sample (although the data points prove 
marginally higher than the model). 

The luminosity function highlights the difficulty in identify- 
ing the low-luminosity systems in the sample. Although our ab- 
solute limit for detection is » lO'^^ergs/s in both fields (Fig.l20]i 
at redshifts of 1.1 < z < 2.0, the identification of such low lumi- 
nosity groups proves difficult at Lx < 4 x lO^^'ergs/s. Including 
tentative detections appears to correct for this, however this 
brings greater risks of including objects falsely identified as clus- 
ters. Additionally, at fainter luminosities/fluxes the volume cor- 
rection (oc l/Q(/x)) increases significantly making our cluster 
counts in the lowest luminosity bin in Fig. l22l more uncertain. 

Thus a disagreement between dn/dz of clusters and the cos- 
mological prediction comes primarily from lowest mass clusters 
or groups close to the detection limit and not from the massive 
clusters, where detection is clear and identification is more obvi- 
ous. A similar picture has been seen in Finoguenov et al. (20ly) 
with the conclusion that, for a robust assesment of the cosmol- 
ogy, one needs to increase the area of such surveys to at least 
50deg^ and use the most massive clusters in the sample. We 
note, that in the SXDF field no such massive clusters have been 
found Finoguenov et al. (2010). It is also quite conceivable that 
the fraction of valid flag=5 identifications increases as the mass 
of the system decreases, which can be verified through the spec- 
troscopic follow-up. 



6. Conclusions 

Using a combination of XMM-Newton X-ray data and opti- 
cal/NIR photometry, we have performed a survey of z > 1.1 
cluster candidates in the CFHTLS Dl and D4 fields. Crucially, 
NIR data is required in order to apply the red-sequence analysis 
at such redshifts (as the 4000A break moves out of the optical 
into the NIR bands). As such the WIRDS data, which provides 
the best combination of wide field and deep NIR data currently 
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Fig. 22. X-ray luminosity function calculated for 1.1 < z < 2.0 
candidate clusters in the CFHTLS Dl and D4 fields. The filled 
circles show the X-ray luminosity function for flag<3 cluster 
candidates only, whilst the open circles show the result for all 
candidates with flag<5. The solid line shows the predicted lu- 
minosity function based on our survey geometry and WMAP 7- 
year cosmological parameters. The dashed line shows the same 
but with the WMAP 7-year Q.^ and cr^ parameters reduced by 
20-. 



available, is one of the best currently available resources for the 
identification of high redshift clusters. Using the WIRDS data in 
conjunction with CFHTLS deep optical photometry, we identify 
atotalof 15z > 1.1 cluster candidates in a total area of l.Odeg. 
All candidates were selected based on X-ray emission in deep 
XMM data, which provides a relatively unbiased sample. 

Of our 15 cluster candidates, 7 are considered firm candi- 
dates (flag<3). These have both strong and clearly extended X- 
ray emission combined with strong red-sequence detections of 
cluster members that are clearly clustered within the X-ray de- 
tection. For a number of these, spectroscopic redshift data is 
available. One of these (WIRDSCS-4-25) is the previously iden- 
tified cluster, XMMXCS J2215.9-1738, at z = 1.45. Our red- 
sequence analysis independently determined a cluster redshift 
estimate of z = 1 .39, showing the success of the red-sequence 
method. We estimate a luminosity and mass for this cluster of 
2.1 X 10'''*ergs/s and 1.9 x lO'^^M©, making it the most mas- 
sive cluster in our sample. Additionally, WIRDS CS-1-76 has 
been previously identified bv lAndreon et al.l (l2009l) . as a z = 1 .9 
cluster (JKS041). As they note, there appear to be a number 
of structures in the line of sight of this X-ray source. Based 
on red-sequence analysis and photometric redshifts we identify 
four apparent structures with redshifts of z = 0.80, z = 0.96, 
Z = 1.13 and z = 1.49. The first three of these are confirmed 
by VVDS Deep spectroscopic data, although the z = 0.80 and 
z - 0.96 structures appear, based on both the spectroscopic and 
red-sequence data, to be dominated by blue star-forming galax- 
ies. The z = 1.13 and z = 1.49 both show strong signs of a clus- 
tered red-sequence correlated with the extended X-ray flux and 
the spectroscopic redshifts for the former again prove the effec- 
tiveness of our red-sequence analysis. Based on our analysis, we 
see no evidence for a cluster at z = 1 .9. Of the remaining flag<3 



cluster candidates, BLOX J2215. 9-1751. 6 has been previously 
identified as a z = 1.1 cluster based on optical photometry only 
(Dietrich et al. 2007; Olsen et al. 2007). With the addition of the 
WIRDS NIR data, we improve upon this redshift estimate, plac- 
ing the cluster at z = 1.17. 

Adding to the 7 high-confidence candidates, we present 8 
further candidates attributed with confidence flag=5, bringing 
our total number of candidates up to 15. These candidates have 
only either weak extended X-ray emission, point-source like X- 
ray detections or non-secure redshift estimates from the red- 
sequence analysis (i.e. due to small numbers of member galaxies 
or poorly clustered member galaxies). By including these candi- 
dates, we place upper limits on the numbers of clusters we are 
able to detect in the given area using the X-ray detection method 
on this data-set. 

Compari ng to other surveys inco rporating z > 1-1.1 clus- 
ter samples, lEisenhardt et al.l (l2008l) reported 106 z > 1 clus- 
ter candidates from the IRAC Shallow Survey (based on optical 
and infrared photometry without any requirem ent for X-ray de- 
tections) with a sky density of x 14.6deg"^. Finogueno v et al.l 
(2007) presented a collection of 8 z > 1 cluster candida tes with a 
sky density of x 3.8deg"^ and lFinoguenov et al.l(l2010h reported 
on 13 z > 1 cluster candidates in the SXDF with a sky density of 
w lO.Odeg"^. Our combination of deep optical and near infrared 
photometry with the XMM X-ray data is therefore an important 
addition to the available surveys of high redshift cluster candi- 
dates. 

We show cluster number counts based on our high redshift 
cluster candidates and compare these to predictions based on the 
WMAP 7-year cosmological parameter estimates. The flag<3 
sample shows lower number counts compared to the WMAP 
7-year model, which takes into account the survey geometry 
and detection method. This tentatively favours lower values for 
the cosmological parameters 0,„ and cr^ than prescribed by the 
WMAP 7-year results and we show that this is comparable to 
a 2cr reduction in both of these parameters. This sensitivity of 
cluster number counts to £!,„ and erg, which is particularly strong 
at z > 1, illustrates the promise of using cluster number counts 
to constrain cosmology. Based on our survey, we highlight the 
issues remaining in applying this method however. As stated, 
modeling the number counts relies on a clear knowledge of the 
detection limits of the cluster survey. We therefore provide a fo- 
cus on non-secure (flag=5) candidates found using our detec- 
tion methods. By including these objects, we find good agree- 
ment between the WMAP 7-year model and the cluster number 
counts. Significantly, we also present the z > 1.1 cluster lumi- 
nosity function for our sample and find good agreement between 
our data and the WMAP 7-year model, except at faint luminosi- 
ties, where we only find good agreement if we include the flag=5 
candidates. We conclude that such non-secure candidates will ul- 
timately be a combination of correct and false cluster detections, 
which we find introduce an uncertainty in our survey equiva- 
lent to the 2cr constraints on Q„, and cr^ from WMAP. Any con- 
straints on cosmology also rely on the cluster scaling relations, 
which are required to constrain the X-ray luminosity limit of 
the survey. We discuss the eff'ect of the relations, noting that the 
scatter on the scaling relations is a key issue and in addition it 
remains for these relations to be calibrated beyond z > 1 . 

The observations presented here suggest that there are per- 
haps too few clusters at z > 1, based on the numbers of rel- 
atively secure identifications of groups and clusters, compared 
to predictions using WMAP-7 cosmology. However, the current 
level of statistical and systematic uncertainties prevent us draw- 
ing a secure conclusion. Ultimately, cluster number counts can 
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provide a complimentary and independent method for constrain- 
ing the cosmological model. In the longer term, the techniques 
presented here (combining X-ray cluster detection followed by 
red-sequence identification) present an unparallelled technique 
to deriving relatively unbiased group and cluster samples. With 
such samples we will be able to independently constrain Q„,, 
cr^ and determine the level of non-gaussianity in the primordial 
density fluctuations. 
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